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A seasonal variation of glucose homeostasis in humans has been reported in various 
geographic regions. In this study, we examined seasonal variations in hemoglobin Ale 
(HbAlc) in patients with type 2 diabetes living in Korea. We analyzed 57,970 HbA1c values 
from 4,191 patients and the association of these values with ambient temperature for 3.5 
yr. Overall, HbAlc exhibited its highest values from February to March and its lowest 
values from September to October (coefficient for cos t = -0.0743, P = 0.058) and the 
difference between the peak and nadir in a year was 0.1 6%-0. 25%. A statistically 
significant seasonal variation was observed in the patients who were taking oral anti- 
diabetic drugs (OADs) without insulin treatment (coefficient for cos t= -0.0949, 
P< 0.05). The Spearman correlation coefficient between daily HbAlc values and the 
corresponding 3-month moving average ambient temperature was -0.2154 (95% 
confidence interval [CI]: -0.2711 , -0.1580; P< 0.05). In conclusion, HbAlc values 
exhibited a seasonal variation in Korean patients with type 2 diabetes, with the highest 
values during the cold season, particularly in those who were treated with OADs, which 
should be taken into account in clinical practice for stable glucose control during the cold 
season. 
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INTRODUCTION 

Because the life of all animals, including humans, is constantly 
influenced by the environment, time-dependent variations in 
biological functions, such as circadian or circannual rhythm, 
are commonly found in animal physiology (1-3). Interestingly, 
seasonal variations are frequently observed in metabolism of 
mammals. The most striking example is the seasonal variation 
in glucose metabolism in the desert gerbil Psammomys obesus 
(sand rat), exhibiting nutritionally induced insulin resistance 
and hyperglycemia (4), which is observed in spring and autumn 
when the animal is in its natural habitat (5). 

Glycemic variation has also been observed in humans. He- 
moglobin Ale (HbAlc), which is the product of non-enzymatic 
glycation of the hemoglobin molecule and reflects the average 
plasma glucose concentrations over the previous 2-3 months 
(6), has been shown to exhibit seasonal variations in both type 



1 and 2 diabetic patients. For example, seasonal variation in 
HbAlc was shown in young type 1 diabetic patients in Poland, 
with the highest values in winter and lowest values in summer 
(7); in children with type 1 diabetes in the U.K., with lower lev- 
els of HbAlc during the summer months (8); in type 2 diabetic 
veterans in the USA, with its peak during March to April and 
trough during September to October (9); and in Japanese dia- 
betic patients with type 1 or type 2 diabetes, with the highest 
values in March and lowest values in August (10). These studies 
were performed mainly in the northern hemisphere. Intrigu- 
ingly, a study encompassing multiple geographic regions, in- 
cluding Calgary, Edmonton, Wisconsin, Singapore, and Mel- 
bourne, revealed that HbAlc levels in subjects with or without 
diabetes were lowest during warmer seasons and highest dur- 
ing cooler seasons regardless of hemisphere (11), indicating 
that temperature or other seasonal factors may influence plas- 
ma glucose levels. Because Korea has four distinct seasons, we 
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hypothesized that ambient temperature or season may affect 
the glycemic control of patients with type 2 diabetes. Therefore, 
in this study, we examined the association between HbAlc and 
calendar day and/or ambient temperature in Korean patients 
with type 2 diabetes who received various types of anti-diabetic 
treatments. 

MATERIALS AND METHODS 
Subjects and data collection 

We retrieved data from electronic medical records of the Seoul 
National University Hospital from October 2007 to May 2011. 
Mean daily temperature data in the Seoul area during the study 
period were obtained from the Korea Meteorological Adminis- 
tration. The data collected from electronic medical records in- 
cluded age, sex, diagnosis, HbAlc, prescribed medications and 
residential addresses. In total, we collected 128,284 HbAlc val- 
ues measured using high performance liquid chromatography 
(Variant II Turbo, Bio-Rad, San Francisco, CA, USA) from 14,689 
patients (Fig. 1). We excluded subjects if: 1) age < 30 yr; 2) pa- 
tients had type 1 diabetes; 3) patients were treated with immu- 
nosuppressants, diagnosed with cancer or had been hospital- 
ized during the study period; and 4) HbAlc measurements were 
made fewer than three times per year. According to these crite- 
ria, 10,498 patients and 70,314 HbAlc observations were exclud- 
ed. Therefore, a total of 4,191 patients (2,211 men and 1,980 wo- 
men, 65 ± 10 yr old) and 57,970 HbAlc observations were sub- 
jected to the final analyses. The patients were classified into four 
categories arbitrarily determined according to baseline HbAlc: 



HbAlc < 7% (n = 1,727); 7% < HbAlc < 9% (n = 2,148); 9% < 
HbAlc < 11% (n = 274); and 11% < HbAlc (n = 42). We also di- 
vided the patients into five groups based upon the type of anti- 
diabetic therapy as follows: a group using lifestyle modification 
alone (n = 166), a group using oral anti-diabetic drugs (OADs) 
(n = 3,022), a group using insulin therapy only (n = 172), a group 
using combined therapy of OADs and insulin (n = 186), and a 
group who changed treatment from OADs to add-on insulin 
therapy or insulin- only therapy (n = 645). 

Statistical analysis 

Statistical analyses were carried out using SAS software (SAS 
Institute, Cary, NC, USA). First, the time series analysis was used 
to find seasonality in HbAlc. Next, we performed a Spearman 
correlation analysis to evaluate the correlation between HbAlc 
and ambient temperature because neither variable followed a 
normal distribution. In the time series analysis, monthly mean 
HbAlc values were determined by averaging the middle 95% of 
HbAlc values over a one month interval after excluding the up- 
per and lower 2.5% of the values. The time series analysis of mon- 
thly mean HbAlc values was modeled as a linear combination 
of time t, cos t, sin t, Ei, and e 2 . The coefficient of t represents a 
linear trend in HbAlc values corresponding to time change. 
The coefficients of cos t and sin t represent a cyclic patterns in 
HbAlc variation, such as seasonality. The noise terms Ei and 
e 2 were derived separately from the first and second autore- 
gressive parts of the time series model and represent the ran- 
domness of the HbAlc values. The total time series model was 
fitted to monthly mean HbAlc and coefficients, and correspon- 



The provisional population (Oct. 2007-May 2011) 
(n = 14,689) 



Excluded patients (n = 10,498) 

1) Patients younger than 30 yr old (n = 351) 

2) Type 1 diabetic patients (n = 56) 

3) Patients who were treated with immunosuppressants (n = 937) or 
diagnosed with cancer (n = 2,016) or had been hospitalized during 
the study period (n = 595) 

4) HbA1 c measurements fewer than 3 times per year (n = 6,543) 



The final population (Oct. 2007-May 201 1) 
(n = 4,191) 



Stratification by the type of anti-diabetic treatment 

1) Lifestyle modification alone (n = 166) 

2) OADs only (n = 3,022) 

3) Insulin therapy only (n = 172) 

4) Combined therapy of OADs and insulin (n = 186) 

5) Changing treatment from OADs to insulin therapy* (n = 645) 



Stratification by baseline HbAlc value 

1) HbAlc < 7% (n = 1 ,727) 

2) 7% < HbAlc < 9% (n = 2,148) 

3) 9% < HbA1 c < 1 1 % (n = 274) 

4) 11% < HbAlc (n = 42) 



Fig. 1. Flow chart of patient selection and grouping, n, number of patients; OADs, oral anti-diabetic drugs. *Add-on insulin therapy or insulin-only therapy. 
http://dx.doi.0rg/10.3346/jkms.2014.29.4.550 http://jkms.org 
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ding P values were determined. Time series analyses were con- 
ducted separately in the whole patient group and subgroups 
divided by types of anti-diabetic therapy or baseline HbAlc val- 
ue. Spearman correlation analysis was performed using a daily 
mean of the middle 95% of HbAlc values as the response vari- 
able. AP < 0.05 was considered statistically significant. 

Ethics statement 

This study was approved by the institutional review board of the 
Seoul National University Hospital (IRB No. C- 1103- 149-357). 
Informed consent was exempted by the board. 

RESULTS 

Time-series analysis 

The HbAlc values exhibited a decreasing secular trend over the 
study period (the coefficient for t = -0.0042, P < 0.05) (Fig. 2A 
and Table 1). Overall, a trend of seasonal variation in the month- 
ly average of HbAlc was observed with highest values from Fe- 
bruary to March and lowest values from September to October 



(the coefficient for cos t = -0.0743, P = 0.058). The differences 
between the highest and lowest HbAlc values (i.e., amplitudes) 
in a year were 0.16%-0.25%. However, a distinct seasonal varia- 
tion in HbAlc was noticed in the subgroup using OADs only 
(the coefficient for cos t = -0.0949, P < 0.05), indicating that gly- 
ceric control was better in the warmer season and worse in 
the cooler season for this group of patients (Fig. 2B and Table 1). 

Correlation between HbAlc and temperature 

The Spearman correlation coefficient between daily mean HbAlc 
values and 3 month-moving averages of daily ambient temper- 
ature was measured (Table 2). Overall, the correlation coeffi- 
cient between these two variables was -0.2154 (95% confidence 
interval [CI] -0.2711, -0.1580; P < 0.05), suggesting that a lower 
ambient temperature may have a negative influence on glyce- 
mic control in patients with type 2 diabetes. In subgroup analy- 
sis, however, this negative correlation was not statistically sig- 
nificant in subjects who were treated with lifestyle modification 
alone or insulin or whose baseline HbAlc was > 1 1% (Table 2). 



Observed mean HbA1 c 
Predicted mean HbAlc 




Observed mean HbAlc 
Predicted mean HbAlc 




^ ^ ^ ^ q 



Fig. 2. Temporal variation in monthly mean HbAlc, observed levels and predicted levels using the time series model in the overall population (A) and in the group taking oral 
anti-diabetic drugs only (B). 



Table 1. The time series analysis in the overall population and subgroups 



Subjects 








sin f § 


cos f § 












e 2 " 


Overall (n = 4,191) 


-0.0042 


± 


0.0010* 


-0.0071 


+ 0.0507 


-0.0743 + 


0.0379 f 


-0.1323 


+ 


0 


1576 


-0.2948 


+ 0.1633 


Subgroups according to the type of anti-diabetic treatment 




























Lifestyle modification alone (n = 1 66) 


0.0006 


± 


0.0020 


-0.0041 


+ 0.1188 


-0.0455 + 


0.0885 


-0.1052 


+ 


0 


1625 


-0.2012 


+ 0.1699 


OADs only (n = 3,022) 


0.0037 


± 


0.0010* 


0.0216 


+ 0.0506 


-0.0949 + 


0.0379* 


-0.1531 


+ 


0 


1569 


-0.2896 


+ 0.1646 


Insulin therapy only (n = 172) 


-0.0025 


± 


0.0011* 


-0.0879 


+ 0.1051 


-0.0022 + 


0.0775 


0.0480 


+ 


0 


1695 


0.0283 


+ 0.1784 


Combined OADs and insulin (n = 186) 


-0.0039 


± 


0.0010* 


0.0592 


+ 0.0823 


-0.1150 + 


0.0607 f 


0.2089 


+ 


0 


1677 


-0.1433 


+ 0.1643 


Changed from OADs to insulin therapy 1 (n = 645) 


-0.0084 


± 


0.0011* 


-0.0910 


+ 0.0877 


-0.0245 + 


0.0647 


0.0024 


+ 


0 


1659 


-0.0291 


+ 0.1702 


Subgroups according to the baseline HbAlc 




























HbAlc < 7% (n = 1,727) 


0.0051 


± 


0.0034 


-0.0444 


+ 0.0547 


-0.0480 + 


0.0418 


-0.4725 


+ 


0 


1534* 


-0.3772 


+ 0.1626* 


7% < HbAlc < 9% (n = 2,148) 


-0.0094 


± 


0.0011* 


0.0023 


+ 0.0553 


-0.0725 + 


0.0415 


-0.2818 


+ 


0 


1643 


-0.1426 


+ 0.1658 


9% < HbA1 c < 11 % (n = 274) 


-0.0305 


± 


0.0138* 


0.0316 


+ 0.1529 


-0.0510 + 


0.1176 


-0.6927 


+ 


0 


1618* 


-0.2171 


+ 0.1678 


11% < HbAlc (n = 42) 


-0.0100 


± 


0.0044* 


0.2970 


+ 0.3925 


0.0338 + 


0.1810 


0.1803 


+ 


0 


1799 


-0.0136 


+ 0.1787 



The data are presented as the mean + SD. *P< 0.05, + 0.05 < P< 0.06. the first autoregressive error; e 2 , the second autoregressive error; OADs, oral anti-diabetic drugs. 
*The coefficient of t represents a linear trend in HbAlc values corresponding to time change. § The coefficients of cos f and sin t represent a cyclic pattern in HbAlc variability. 
"The noise terms ei and e 2 represent the randomness of HbA1 c values. ^Add-on insulin therapy or insulin-only therapy. 
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Table 2. Spearman correlation coefficient between 3-month moving averages of daily 
mean HbA1 c values and daily ambient temperature 



Subgroups 



Spearman correlation coefficient 
(95% CI) 



Overall 


-0.2154 


(-0.2711, 


-0.1580)* 


Subgroups according to the type of anti-dia- 








betic treatment 








Lifestyle modification alone 


-0.0306 


(-0.0965, 


0.0356) 


OADs only 


-0.2122 


(-0.2680, 


-0.1548)* 


Insulin therapy only 


-0.0177 


(-0.0825, 


0.0473) 


Combined OADs and insulin 


-0.0702 


(-0.1333, 


-0.0064)* 


Changed from OADs to insulin therapy + 


-0.1320 


(-0.1897, 


-0.0732)* 


Subgroups according to the baseline HbAlc 








HbAlc < 7% 


-0.1863 


(-0.2428, 


-0.1284)* 


7% < HbA1 c < 9% 


-0.1469 


(-0.2043, 


-0.0883)* 


9% < HbAlc < 11% 


-0.0663 


(-0.1272, 


-0.0048)* 


11% < HbAlc 


-0.0416 


(-0.1474, 


0.0653) 



*P< 0.05; f Add-on insulin therapy or insulin-only therapy. CI, confidence interval; 
OADs, oral anti-diabetic drugs. 



DISCUSSION 

We observed a trend in seasonal variation in monthly averages 
of HbAlc with the highest values in February- March (late win- 
ter-early spring) and the lowest values in September-October 
(early autumn) in Korean patients with type 2 diabetes. Because 
the HbAlc value mirrors the average plasma glucose levels dur- 
ing the previous 2-3 months (6), these peak and trough HbAlc 
values reflect the mean plasma glucose levels during the winter 
and summer seasons, respectively. These results are consistent 
with previous reports showing maximum values of HbAlc in 
cooler seasons and minimum values in warmer seasons (7-12). 
In addition, the differences between the highest and lowest Hb- 
Alc values (i.e., amplitude) in a year were 0.16%-0.25%, which 
is also similar to the results of other studies, which reported 
amplitudes of 0.13%-0.45% (7, 9). 

We examined the interaction between seasonal variation of 
HbAlc and types of anti-diabetic therapy or baseline HbAlc 
levels. Seasonal variations in HbAlc were observed only in pa- 
tients who were treated with OADs, and this was not affected by 
baseline HbAlc values. It is of interest that HbAlc was not af- 
fected by ambient temperature in patients who were treated 
with either lifestyle modification alone or insulin. Perhaps this 
lack of effect means that glycemic control in relatively well- con- 
trolled type 2 diabetes is not affected by ambient temperature. 
In addition, patients using insulin therapy may effectively re- 
spond to the altered glycemic control by adjusting their insulin 
dosage. However, seasonal variations in HbAlc were also re- 
ported in healthy women (13) or patients who were treated with 
insulin (7, 9). Therefore, we cannot guarantee that patients who 
are treated with either lifestyle modification alone or insulin are 
less susceptible to seasonal changes in glycemic control. Fur- 
thermore, because the number of the patients in these subgroups 
was relatively small, we may have missed the actual signal of 



seasonal variation. 

Although seasonal variations in HbAlc have been reported 
by several independent researchers (7-11), the mechanism un- 
derlying this phenomenon is still not definite. Our study was 
primarily concerned with the correlation between HbAlc and 
ambient temperature and the impact of sun exposure or differ- 
ent levels of physical activity in glucose homeostasis has not 
been considered. Nonetheless, serum 25-hydroxyvitamin D 
concentration, which is affected by exposure to ultraviolet light 
(14), or serum melatonin concentration, which is dependent 
on photoperiod (15), could play an important role. It was re- 
ported that vitamin D deficiency may increase insulin resistance 
and the risk of metabolic syndrome (16-18). Therefore, a lower 
serum level of 25-hydroxyvitamin D during the winter season 
(19) may aggravate glycemic control. Melatonin production is 
highest at night (20), and some studies reported seasonal varia- 
tions in human melatonin production with increased levels in 
winter (15, 21). It was reported that the nighttime increase in 
melatonin was associated with rises in circulating glucose lev- 
els and reduced insulin sensitivity (22). However, there are stu- 
dies showing that a lower level of melatonin was associated with 
development of type 2 diabetes (21, 23). Therefore, further stu- 
dies are needed to examine the role of melatonin in the season- 
al variations in HbAlc particularly in patients with type 2 dia- 
betes. 

Aside from seasonal biological variations, cultural events or 
food availability may influence the seasonal variability in HbAlc. 
In a Chinese study examining HbAlc during the Chinese New 
Year's holiday, mean HbAlc was increased by 0.094% ± 0.828%, 
most likely due to decreased physical activity, increased intake 
of foods and/or increased consumption of alcoholic beverages 
(24). In addition, there is a tendency for increased intake of ani- 
mal fat-containing foods and decreased physical activity in win- 
ter (25-27). Furthermore, it was reported that body mass index 
and waist circumference were higher in winter than in summer 
(28), which could be explained by changes in calorie intake and 
physical activity. These non-biological factors should be con- 
sidered when evaluating seasonal variations in HbAlc because 
these factors could be modified by education and counseling. 

Of interest, there was a decreasing secular trend of monthly 
averaged HbAlc values over time during the study period, as 
depicted in Fig. 2. This trend reflects the movement towards 
stricter glycemic control for patients with type 2 diabetes, which 
was supported by compelling evidence for the "legacy effect" of 
intensive glycemic control from the initial diagnosis of type 2 
diabetes based on long-term follow-up in the U.K. Prospective 
Diabetes Study (29). Additionally, the Korean Diabetes Associa- 
tion recommended a target HbAlc level of < 6.5% in Novem- 
ber 2007 (30), and this more stringent target level may have re- 
sulted in an improvement in HbAlc values during the study 
period. 
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There are some limitations and strengths in our study. As men- 
tioned earlier, we did not measure the level of sun exposure or 
physical activity, which may affect glucose homeostasis in pa- 
tients with type 2 diabetes. The lack of data on 25-hydroxyvita- 
min D and melatonin prohibited us from examining the associ- 
ation between these hormones and glycemic control. However, 
the strength of our study is that we included more than 4,000 
type 2 diabetic patients with HbAlc values measured on a regu- 
lar basis by a standardized method in a single institution over 4 
yr, which allowed us to examine the seasonal variation accord- 
ing to the different types of anti-diabetic therapy. 

In conclusion, we observed a trend in seasonal variation in 
HbAlc in Korean patients with type 2 diabetes, particularly in 
the OAD-treated patients. Although the mechanism underlying 
this phenomenon is yet to be elucidated, stricter adherence to 
lifestyle modification and frequent assessment of glucose con- 
trol appear to be necessary during the winter season to prevent 
aggravation of hyperglycemia in patients with type 2 diabetes. 
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